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Abstract. We study Jjij/ production in pp collisions at a/s = 1.96 and 7 TeV using 
the Colour-Singlet Model (CSM), including next-to-leading order (NLO) corrections and 
f~^ dominant a^ contributions (NNLO*). We find that the CSM reproduces the existing data if 

CN the upper range of the NNLO* is near the actual -but presently unknown - NNLO. The direct 

r^ yield polarisation for the NLO and NNLO* is increasingly longitudinal in the helicity frame 

>_Z when Pj gets larger Contrary to what is sometimes claimed in the literature, the prompt 

I Jlijj yield polarisation in the CSM is compatible with the experimental data from the CDF 

collaboration, when one combines the direct yield with a data-driven range for the polarisation 
'^ ofJIipfmmxc- 

p , 

53^ 1. Introduction 

The numerous quarkonium-production puzzles at hadron colliders were attributed not too 
^ long ago to non-perturbative effects associated with channels in which the heavy quark- 

C^l pair is produced in a colour-octet state [1]. a^ and or^ corrections to the CSM [2] are now 

psj widely recognised as essential to understand the Pj spectrum of J I if/ and T produced in high- 

"^ energy hadron collisions [3, 4, 5, 6, 7]. This calls for a factorised description of high-Pr J I if/ 

Pr beyond leading power [8]. The effect of QCD corrections is also manifest in the polarisation 

^H predictions. While the J/ifj and T produced inclusively or in association with a photon are 

predicted to be transversely polarised at LO, it has been found that their polarisation at NLO 
is increasingly longitudinal when Pj gets larger [5, 6, 9]. In recent works [10, 11], we have 
also shown that the CSM alone is sufficient to account for the magnitude of dcr/dy at RHIC, 
^ Tevatron and LHC energies. 

We evaluate here the Pj dependence of the J /iff yield and its polarisation at Tevatron and 
LHC energies. We describe the procedure used to obtain a first evaluation of some dominant 
contributions at a^ (NNLO*) in addition to the yield at NLO (up to a^). We then compare 
available data from the Tevatron and the LHC with our results: the direct yields differential 
in Pt along with the polarisation vs Pt for the prompt yield using an essentially data-driven 
estimation of the polarisation for J/ifr from;^,,- 



X 



J I ifj production at yfs = 1.96 andl TeV 2 

2. Cross-section 

For the NLO cross section, we use the partonic matrix elements of [3]. In order to investigate 
the expected impact of NNLO QCD corrections for increasing Pj, we also present the NLO 
results plus the real-emission contributions at a^ evaluated along the lines of [6], referred 
to as NNLO*. At a^, the last:]: kinematically-enhanced topologies open up, with a P^'^ fall 
off of dcrldP\. The procedure used here for the NNLO* is exactly that of [6]: the real- 
emission contributions at a^ are evaluated using MADONIA [12] by imposing a lower bound 
on the invariant-mass squared of any light partons {Sij). The dependence on this cut should 
decrease for larger Pj since no coUinear or soft divergences can appear there for the new 
channels opening up at a^ with a leading-Pj behaviour, i.e. the ones which interest us. For 
other channels, whose Born contribution is at or^ or or^, the cut would produce logarithms 
of Sjj/s^^". These are not necessarily small, but they are expected to be factorised over their 
corresponding Bom contribution, which scales as P~^ or P~^. They are thus suppressed by at 
least two powers of Pt with respect of the leading-PT- contributions (Pj-'^). The sensitivity on 
s'^'^ is expected to be small at large Pj- 
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Figure 1 : daldPr x Br for direct Jjifj production from NLO and NNLO* CS contributions at 
V^ = 1.96 TeV (left) and dX ^fs = 1 TeV for central (middle) and forward (right) rapidities. 
These are compared to the CDF [14], ATLAS, CMS and LHCb data [15, 16, 17] multiplied 
by a constant direct fraction from CDF [18]. See text for details on theoretical-error bands. 



Our results are shown on Fig. 1. The CSM is very close to the existing data, if 
the upper range of the NNLO* is a relevant evaluation of the NNLO. The uncertainty 
bands at NLO are obtained from the combined variations of the charm-quark mass (m^ = 
1.5 + 0.1 GeV), the factorisation Hf and the renormalisation hr scales chosen in the couples 
((0.75, 0.75); (1, 1); (1, 2); (2, 1); (2, 2)) x m^ with m^ = 4m^ + P].. The band for the NNLO* 
is obtained using a combined variation of m-c, O.Smr < /ur = fip < 2mT and 2.25 < i'J"'" < 9.00 
GeV^. We have used the NLO set CTEQ6_M [13] and have taken \Rj/^{0)\^ = 1.01 GeV^ and 
Br(7/i/r ^ rr) = 0.0594. 

i We do not expect any further kinematical enhancement as regards the Pt dependence when going further 
in the as expansion: P^* is the slowest possible fall-off. Above a^ , usual expectations for the impact of QCD 
corrections would then hold. One would expect a K factor multiplying the yield at NNLO accuracy, which would 
be independent of Pj- and of a similar size as those of other QCD processes. A farther enhancement by an order 
of magnitude between the NNLO and N^LO results would be quite worrisome. 
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3. Polarisation 



The polarisation parameter a is extracted bin by bin in y or Pt from the normalised distribution 
of the polar angle 9 between the ^^ direction in the J/tff rest frame and its direction in the 
laboratory frame, /(cos 6) = ^, ^^3, (1 + a cos^ 6). We thus work in the helicity frame, a is also 
related to a ratio of the polarised cross sections: a = ""^'l"''- . 

For the time being, there does not exist any measurement of direct J/i/r polarisation, 
i.e. after the extraction of the^^ feed-downs (up to 30-40 %) which may strongly impact on 
the observed values of a. It is however possible to constrain its effects by using existing data 
on cr^^Jcr^^^ and by relying on El dominance for the transition ;^c. -^ J I if/ -\- y. 

Indeed, using El dominance [19], one can obtain [21] a range of the yield of 
longitudinally (transversely) polarised Jjif/ in terms of simple relations involving the polarised 
Xc yields. Allowing for extreme cases, these relations allow the yield from Xc to be fully 
transversely polarised, while there is a minimal value of a. Following the discussion of [21] 
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Figure 2: Comparison between the extrapolation of a for prompt J I if/ in pp at ^fs = 1 .96 TeV 
(red band), the direct NLO a (gray line), the direct NNLO* a (thinner dark-red band) and the 
CDF data [22] for prompt J/ if/. 

Since 30% of the J /iff come from Xc nearly independent of Pt in the range considered 
here [18], we expect a partial contribution to the polarisation ranging from 0.3 x (-1-1) to 
0.3 X (-0.42). Regarding the other 70%, one multiplies the result for the direct yields by 
0.7, since the polarisation of J/if/ from if/(2S) is expected to be identical to the direct one. 
Doing so, one obtains the extrapolation shown on Fig. 2. If the Jfifr from;^^^ yield is strongly 
transversely polarised (the upper limit), the polarisation of the prompt yield is in rather good 
agreement with the data. 

4. Conclusion 



We have evaluated the NLO and NNLO* J /iff yield at Tevatron and LHC energies. As found 
for T at the Tevatron [6] and for Jfif/ at RHIC [21], the upper bound of the CSM predictions 
is very close to the experimental data from CDF, ATLAS, CMS and LHCb. However, the 
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NNLO* evaluation is not a complete NNLO calculation. It is affected by logs of an IR cut- 
off whose effect might not vanish as quickly as one has anticipated. It may very well be 
that the upper limit of the prediction -close to the data- accurately reproduces the complete 
NNLO yield, or that the lower limit -close to the NLO yield- reproduces the NNLO yield. If 
the upper limit of the NNLO* does indeed overestimate the NNLO, the CSM alone is likely 
insufficient to account for the data. Conversely, the CSM alone is enough and the colour-octet 
contributions are not required. 

As regards polarisation, we have derived a range for the prompt yield polarisation. 
This range is affected by admittedly large theoretical uncertainties, but the upper edge - 
corresponding to a transversely polarised feed-down- is in rather good agreement with the 
data from the Tevatron. We recall that the trend for a longitudinally (direct) i/r(25 ) yield [22] 
was also met by the NNLO* [7]. 

In conclusion, the CSM may very well provide a good description of J/iff production in 
high-energy pp collisions, both in terms of the cross section and the yield polarisation. 
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